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Summary. — The two-dimensional relativistic configurational r-space is proposed 
and the exactly solvable finite-difference model of the harmonic oscillator in this 
space is constructed. The wave functions of the stationary states and the corre- 
sponding energy spectrum are found for the model under consideration. It is shown 
that, they have correct non-relativistic limit. Explicit form of the plane momentum 
operator is determined. 
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1. Introduction 

The harmonic oscillator is one of the important problems of the non-relativistic quan- 
tum physics. It is one of few problems of the quantum mechanics having an exact 
solutions. The theory of the harmonic oscillator has various applications in the atomic 
and molecular as well as nuclear and particle physics [1]. 

Usually, analytical solutions of the one- and three-dimensional non-relativistic quan- 
tum oscillator models are well known due to their enormous importance for expla- 
nation of the various quantum physics phenomena. The solution of the stationary 
Schrodinger equation for the one-dimensional non-relativistic oscillator is simple leading 
to the wave functions expressed through the Hermitc polynomials, whereas the similar 
three-dimensional solution leads to the radial wave function that is expressed through the 
generalized Laguerre polynomials. One need to note here the two-dimensional model of 
the non-relativistic quantum harmonic oscillator, that also has an explicit analytical so- 
lution for the stationary states similar to the three-dimensional model [2] and this model 
has a lot of important applications, too. For example, in [3], the two-dimensional model 
of the harmonic oscillator is used to study the Aharonov-Bohm effect. This oscillator 
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model is also applied to describe a nonlinear amplifier, exemplified by a two-junction 
superconducting quantum interference device in the presence of thermal noise [4] and 
its perturbed model is employed to obtain the Rayleigh-Schrodinger energy perturba- 
tion series [5]. One can note a recent interesting work [6], that considers different two- 
dimensional oscillators, namely a nonlinear Lachmanan-Mathews oscillator, a nonlinear 
oscillator related to the Riccati equation and the standard harmonic oscillator on con- 
stant curvature spaces as nonlinear deformations of the linear harmonic oscillator as well 
as another quite interesting work [7] uses a variational method with a normal ordering 
technique to express the Hamiltonian of two-dimensional anharmonic oscillators, in terms 
of Wicks ordered-product of the functions of the generalized creation and annihilation 
operators in order to find the eigenvalues of the 2D anharmonic Hamiltonian. 

However, almost all of these studies has been done in the framework of the non- 
relativistic quantum approach. To consider mentioned above the two-dimensional prob- 
lem in the framework of the relativistic quantum approach, one of the ways is to solve 
Dirac or Klein-Gordon equations analytically. An approach to solve the Dirac equation 
for quantum harmonic oscillator, that proposed in [8] allows to consider Dirac oscillator 
in two-dimensions, too [9], where the problem has explicit solutions under requirement 
that the potential should be introduced as a pseudovectorial term shifting the momentum 
operator. However, this is not only the way to construct the two-dimensional relativistic 
oscillator model, and another elegant relativistic approach is the relativistic quantum 
mechanics based on the finite-difference equations. 

The purpose of this paper is to construct the relativistic model of the two-dimensional 
harmonic oscillator [2]. Our construction is based on the relativistic finite-difference 
quantum mechanics, developed in [10-20]. Here, the relativistic configuration space is 
the key conception [14]. Note that, the models of the linear, three-dimensional and 
TV-dimensional harmonic oscillators already have been studied in the framework of the 
relativistic finite-difference quantum mechanics [10-13,22]. 

The structure of the paper is as follows: We consider two-dimensional relativistic 
configurational r (x, y)-space in the second section. The model of the harmonic oscillator 
in this space is proposed in the third section. 

2. Two-dimensional relativistic configurational r (x, y)-space 

Two-dimensional relativistic equation for the wave function in the momentum p(p x , Py)~ 
representation in the case of the equal masses has the following form [11, 14]: 



(1) (E p -E)i>(p) = —2 / V U k; E) $ (k) dQ k , 



where dfl^ = dkj y 1 + k 2 /m 2 c 2 is an invariant volume element of the integration in 

the ^-representation, E p — cy 1 p 2 + m 2 c 2 is the relativistic energy and V \ p, k; Ej is the 
interaction potential. 

We carry out the integration in the over the mass shell of the particle with mass m, 
i.e. over the upper sheet of the hyperboloidp 2 ,— p 1 = m 2 c 2 (two-dimensional Lobachevsky 
space). 

Transition to relativistic two-dimensional configuration r-space 
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(2) 



is performed by the use of expansion on the matrix elements of the representations of 
the two-dimensional Lobachevsky space: 



(3) 



Pa - pn 



mc 



-l-ir/X 



where r = rfi, < r < oo, n = (cos ip, sin ip, ) , po — E p /c and X = h/mc is the 
Compton wavelength of the particle. 

As a result of this transition, the eq. (fT]) in the case of the local potentiald of the 
interaction has following finite-difference form: 



(4) [H + V(r)}i,(r) = Eij>(f) 

with the free Hamiltonian Hq of the form: 



(5) 



Hn = mc 



COsh {lXd r ) + ^ Sinh (lXd r ) - ( . 



q2 n iXd r 



Here, the momentum operator of the free particle has following explicit determination: 



(6) 



p = mc ■ n \ — e r 

mc z 



h 



JXd r 



r + iX/2 

that is expressed by emploing the two-dimensional vector m of the following components: 

m = i{— sin ip ■ d v , cos ip ■ d v ) . 

Some useful relations for the operators H, rh and L — ~id v are presented in the Ap- 
pendix. The components of the p momentum vector and free Hamiltonian Ho commutate 
as follows: 



[Pi, Pi J 



Pi, H Q 



0, i,j = x,y. 



The function ([3]) being the relativistic 'plane wave' has the property to be common 
eigenfunction of the operators Hq and p, i.e. 



(7) H Q ■ £(p,r) = E p ■ £(p,r) , p ■ £ (p, f) = p ■ £ (p, r) . 

The relativistic plane waves ([3j) form orthogonal and complete system of the functions: 
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(8) 



{2ttT%) 



"772 / C{i>,r) ■ t(tf,r)w(r)dr = 6 (p (-)&), 
m) J 



where, w (r) = - (— r/A)^ 1 ^ is the weight function [22], 5(p(— )p > ) 
and is the generalized degree [14] having the form: 



(*) = , { r(-ir/A+j) 
A/ ~* r(-ir/A) ' 



(s) 



In the non-relativistic limit we have the following reducing relations: 



lim (Ho - mc 2 ) = H?« = -h. ^ + \ - , 
(9) lim p = Pnr — — ft (indr + m- J , 

c^oo \ 7* y 

lim £ (p, r ) = e 4p " 

3. The relativistic finite-difference model of the two-dimensional harmonic 
oscillator 

We consider the model of the relativistic two-dimensional oscillator that corresponds 
to the following potential of the interaction: 



(10) V{r) = ^f- ^§j- 2 [r (r + iX) XH ■ d%\ e«* . 

It is clear that in the non-relativistic limit we will have: 

lim V (r) = -mbj 2 r 2 . 

c^oo 2 

Due to commutativity of the operators H = H + V (r) and L = —id v , eigenfunctions 
of the Hamiltonian H depend on the angle ip by standard manner: 



(11) 



V(r)= (-r/X) 



(1/2) 



R\m\ (r) 



a tmip 



2tt' 



m = 0,±1,±2, . . 



Therefore, considered two-dimensional problem will be reduced to the linear one, i.e. 
we have to find the eigenfunctions and eigenvalues of the radial part of the Hamiltonian 
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(12) 



H ■ R\ m \ 0) = E\m\ ■ R \m\ (?") 



with the boundary conditions R\ m \ (0) = R\ m \ (oo) = 0. Here H has following explicit 
finite-difference form: 



(13) H = mc 2 { cosh(ad r ) + 



In terms of the dimensionless variable p — r/X and parameter loq — haj/mc 2 , the eq. 
T2|) will have the following form: 



(14) 



1 2/ 

2 W ° I' 



cosh (id p ) + ■f^e* " + ^<4 ( /''"' - ' ) < 



i(i„ 



E\ 



R\m\ (p) = — ^R\m\ (P) , 



with a = m 2 — 1/4 and 7 = bm 2 . 

According to the works [10,12,23], we extract multipliers (— p) and M u (p) — 
lo 1 q P T [ip + v), which determine the behaviour of the radial wave function R\ m \ (p) at 
corresponding points p = and p = 00, and therefore we look for solution of the (|14p as: 

(15) R M ( P ) = C\ m \ ■ {~p) (a) ■ M v (p) • fi (p) , 

where the constants a and v equal to: 



(16) 



(17) 




Then, the function ft (p) will satisfy the following finite- difference equation: 



(18) \(a + ip) (y + ip) e~ l9 » - (a - ip) (v - ip) e i9 »] (p) = 2i-^±VL (a) . 

nu 

Polynomial solutions of the eq. (jXS]) , corresponding to the energy spectrum values 



(19) E\ m \ = E n \ m \ = flu (2n + a + v) , n = 0,1,2,... 

are expressed through the continuous dual Hahn polynomials [23] : 
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(20) n(p) = Sl n (p) =S„(p 2 ;a,M/2). 

The condition (jT9j) gives the quantization rule for the energy level of the two-dimensional 
harmonic oscillator (|10[) and leads to the following form of the radial part of the wave 
function: 



(21) Rn\ m \ (p) = C nM • (~p) {a) ■ M v (p) ■ S n (p 2 -a, is, 1/2) , 

c - / 

™ |m| y n\ ■ r (n + a + v) ■ T (n + a + 1/2) • T (n + v + 1/2) ' 
The functions i? n i m i (p) are orthogonal and normalized by following manner: 



(22) / R n \ m \ (p) ■ R n '\ m \ (p) ■ dp = <W- 



Therefore, from (|22| we obtain following orthonormalization condition for the full 
wave function pip : 



(23) / ^*| m | (f) ■ ip n >\ m >\ (r) ■ w (r) df = 

In the non-relativistic limit, we find that 



1 1 , 1 

hm a = - + -VI + na = -+ m , 

c^oo z z z 



lim [ f 



2' 



(24) 

lim w (r) = 1, 



lim ^S„(e 2 M,;a,^l/2) =L M (^) . 



Here, p 2 = £ 2 /wo, £ = r y /s w i anc ^ are * ne generalized Laguerre polynomials: 

K (x) = ^^- 1 F 1 (-n;p + l,x). 

Consequently, in the non-relativistic limit c — > oo, wave functions (|11[) and energy 
spectrum (|19j) coincide with the corresponding wave functions and energy spectrum of 
the non-relativistic two-dimensional harmonic oscillator, i.e. 
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(25) 



lim </>„ H if) = f) = ■ £ |m| e-^ 2 • LjH (£ 2 ) 



lim (£ nH - mc 2 ) = fiX,, = foj (2n + |m| + 1 



The wave functions ^Qml (^0 are othonormalized as follows: 

(26) | (r)) * • | (0 df = 6 nn , 5 mm , . 

From (|26p we can see that normalization constant has following form: 



C 



NR. 
nlm 



2mujn\ 



1 V hT(n+ \m\ +1)' 



Appendix A. 



The operators n, m and L = — id v satisfy following relations: 



















n x ,L 


272y , 


fly ) ^ 




m Xl L 


= —im y , 


fTly, L 


= im x , 



fii,L' 



(rii + 2i ■ mi) , 



mi,L 2 



mi + 2i ■ mL' 



[fly, m x \ = [n x ,m v l = in x n y , 



[m x ,m y ] = —iL, [m x , n x ] = in yl [%,%] = m 2 



Their scalar products are: 



ft ■ fh = 0, fh ■ n — i, fh = L 2 . 



Under the Hermitian conjugation with respect to the scalar product J (jf ■ ipdr, we 
ve nf — rii, mf — irii — mj, i = x,y. 
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